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Preface 



This report is a collection of material that has been used 
in courses on search, detection and localization modeling. Its 
organization follows to some extent material by S. M. Pollock in 
Selected Methods and Models in Military Operations Research which 
is listed in the report bibliography. The report is not intended 
to be a text on these subjects. In particular, in some areas it 
does not provide the depth of coverage that is found in the book 
Search and Detection by Alan R. Washburn which is cited in this 
report as Reference 21. 
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I. Detection Models and Signal Detection Theory 

Signal detection theory is the basis for analyzing the 
detection models that are described in this report. In signal 
detection theory, the decision making portion of a detection 
system is called the receiver and a detection experiment is the 
observation by a receiver of input data generated during some 
time interval. The data that is related to a target is called 
signal. The data that is not related to the target is called 
noise. In general, the target data is associated with a 
localization region that in some cases is called a resolution 
cell. When a detection experiment is performed, either the event 
H 0 = {the receiver input is noise} or its complement 

= {the receiver input is signal and noise} will occur. In 
the first detection models that are described here, after 
analysis of the input data by a receiver, either the event 
D 0 ={the receiver decides the input is noise} or its complement 
D^ = {the receiver decides the input is signal and noise} will 
occur. Detection models for which D^^ is the complement of D 0 
are called binary detection models or forced choice detection 
models. Four events which are important in binary detection 
model are indicated in the Venn diagram of Figure 1. 

The Venn diagram emphasizes a decision problem that is 
associated with a receiver that can be modeled using a binary 
detection model. The problem is this: Under what conditions 
should the event D^ occur? That is, under what conditions 
should a receiver decide that the input data indicates a target 
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was present in an observed localization region either at the time 
or prior to the time of the observation? 



H 0 H X 

D 0 



Dl 



Figure 1. Four events of importance in binary detection models. 

In the detection model descriptions that follow, the 
following notation and terminology are used: pf = P(Di|Hq), the 

probability of Dj given H 0 , is called the false alarm 
probability; p d = P(D 1 |H 1 ), the probability of given 

is called the detection probability and P = PfHj) , the 
probability of , is called the prior probability. It is the 
probability that a target will be in the localization region when 
the input data is generated. 

In the detection models, the input to a receiver consists of 
a sequence yi/-**ym that is a realization of a discrete 
parameter stochastic process whose random variables are elements 
of a continuous parameter stochastic process. Although the 
receiver input in some cases can be considered to be described by 
a continuous stochastic process, because of the finite amount of 
information (unique data) contained in a bounded sequence of 
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a continuous stochastic process, because of the finite amount of 
information (unique data) contained in a bounded sequence of 
finite length, a discrete stochastic process is sufficient to 
represent the input in these cases. This is established formally 
by the stochastic sampling theorem. In the models, for an 
observation, the input stochastic process has the following 
characteristics: It is a noise process when there is no target 

data and it is a noise process plus a signal process when there 
is target data. In the first detection model that is described 
in Section III, the signal process is a deterministic process. 

In the second and third detection models that are described 
there, the signal process is a random process. To specify the 
noise process or a random signal process, one needs only to 
specify the joint distribution of the finite sequence of random 
variables that determine the process. If the signal process is a 
deterministic process, the signal values are a sequence of values 
that are assumed to be known before an observation is performed. 
Consequently, in this case, to specify the process, one needs 
only to specify these values. 
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II. Decision Criteria 

To simplify the discussion of decision criteria and 
decision rules, a receiver's input will be assumed to be 
determined by a single decision random variable Y. In this 
case, the input process in determined by the conditional 
distribution function Fy(y|Ho) when the input is noise alone 
and by the conditional distribution function FyCylHQ^) when the 
input is signal plus noise. 

The condition that a receiver's input is required to 
satisfy in order that the event D! will occur can be specified 
in terms of a decision rule. For the assumed case, a decision 
rule is a rule which determines for every observable value of Y 
the decision that the receiver is to make. The decision rule can 
be considered to be a function <£(y) which relates each 
observable value y to one or the other of the following two 
commands: d 0 = "Decide that the receiver input was noise." and 

d-L = "Decide the input was signal and noise.". Choosing a 
decision rule 0(y) defines a set n of observable values of Y 
such that the event D^ = { Y e n } . 

The problem which was considered in Section I can now be 
restated in the following way: What criterion should be adopted 

in order to determine a decision rule or, equivalently, its 
corresponding set n ? A desirable characteristic for a criterion 
is suggested by the following argument: Consider the odds in 

favor of given y is observed. That is, consider 
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P(H]Jy = y)/P(H 0 |Y = y) . One might expect that y would be a 
member of the set ft if and only if y made this ratio equal to 
or greater than some value k. But this is equivalent to defining 
n as follows: ft = { y : L(y) > K } where L(y) is the 
likelihood ratio associated with an observed value y and K is 
a constant related to the constant k. This suggests that 
choosing an optimum criterion is equivalent to choosing an 
optimum value for K. Four specific decision criteria are defined 
next in terms of K. For each criterion, n has the above 
form. But for each criterion the choice of K is different. 

The decision criteria are: 

1. The Neyman-Pearson Criterion: Choose ft so that p^ 

is a maximum subject to the constraint that pf < a where a is 
a specified value. For a continuous decision random variable, 
the constant K is chosen so that pf = a. 

2. The Bayes Criterion: Choose fi so that the expected 

cost of a receiver's decision is a minimum. For a continuous 
decision random variable, K = [ (c 10 -c oo)/ ( c 01 ” c ll) Hl“ p )/ p if 
c 10 > c oo an< * C 01 > c ll where c^j is the cost of n Hj . 

3. The Ideal Observer Criterion: Choose ft so that the 

probability that the receiver makes an incorrect decision is a 
minimum. K = (1-P)/P for a continuous decision random variable. 

4. The Minimax Criterion: Choose n when P is unknown so 

that the maximum expected cost of a receiver's decision is a 
minimum. If c^q > c 00 and C 01 > c ll' then 
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K = C ( C 10 -C 00) / ( c 01 -c ll) ]( l -p *) /P* for a continuous decision 
random variable. Here, P* is the value of the prior probability 
P which would make the expected cost of a receiver's decision a 
minimum for this choice of K. 

If a model which specifies the conditional 
distributions Fy(y|H 0 ) and Fy(y|Hi) and a decision rule are 
adopted, then the value of pf and the value of p<j are 
determined. This pair of values (Pf/Pd) is called a receiver 
operating point. If the decision rule results from using a 
likelihood ratio criterion such as one of the four listed above, 
then it will involve the parameter K through the relation 
f] = {y: L(y) > K}. And, for a given value of K, since n 
uniquely determines the pair (Pf /Pd ) • a single operating point 
results. By varying K, a set of operating points can be 
generated which determines a receiver operating characteristic 
curve or ROC curve. Different ROC curves can be produced by 
changing either one or both of the conditional distributions 
which implies changing either the signal process or the noise 
process . 

A decision rule which results from using a likelihood ratio 
criterion in a model in which the input process is determined by 
a set of m random variables can be expressed in terms of a set 
n as follows: n = ( (yi,...,y m ) : L(y!, . . . ,y m ) > K } where K 
is specified in the same way as in the case in which m = 1. 
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III. Three Binary Detection Models 

Three detection models are examined in this section. 

For the first two detection models, the input stochastic process 
for an observation is defined by a time sequence of continuous 
random variables. The random variables represent a sample from a 
continuous parameter stochastic process which is sampled at times 
such that the random variables are independent. For the third 
detection model, the input stochastic process is a counting 
process and it is defined by a single discrete random variable 
that is equal to the number of events that are counted during the 
observation. 

Model I: In the first detection model, a sampled noise 

value is a value of a normally distributed random variable with 
mean zero and with known variance o J and a sampled signal value 
is a known value of a deterministic variable. Thus, the input 
process corresponding to an observation consists of some number 
m of independent normal random variables Y 1 ,...,Y In each with 
variance o’. And, for i = 1,2,3, , m, when a signal is not 

present the mean of Y^ is zero and when a signal is present the 
mean is s^. The result of using a likelihood ratio decision 
rule in the model can be expressed in terms of a random variable 
Z. This random variable is called a crosscorrelation statistic 
and it is defined by Z = E Sj_*Yj. where the sum index 1 = 1,2, 
,m here and in the remainder of this section. However, it 
is more convenient to express the result in terms of a statistic 
V which is defined by V = Z/o 2 . In terms of this random 
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variable, the conditional probabilities pf and p^ are given 

by: pf = 1 - *(v*) and p^ = 1 - 9 (v* - d^) where # 

symbolizes the standard normal cumulative distribution function, 

the threshold v* = (l/o z ) (o’ln K + (1/2) E s{) and 

d = E s i/a*. The parameter d is called the detection index. 

Often, the input stochastic process represents a quantity 

whose square is proportional to power. In such a case, the 

, 2 

average receiver input power is the random variable E Yj/m. If 
a signal is not present, the expected average receiver input 

power is N = E a*/m = a* where N is called the noise. The 

. . . . 2 

average receiver input signal power is S = E s-jyrn where S is 

called the signal. In these terms, d = m* (S/N) where S/N is 

called the signal-to-noise ratio. 

If a receiver's input can be considered to be a time 

sequence of continuous voltage values such as in the case of a 

sonar receiver, in some cases a frequency representation can be 

used that involves the concept of receiver bandwidth. In these 

cases, the noise process is assumed to be such that m = t/St 

where t is the integration time (the duration of an 

observation) and St is the time between samples with 

St = 1/ [ 2 ( BW ) ] where BW is the bandwidth and St is 

determined by the sampling theorem. This implies that the 

detection index can be written as d = 2t* (BW) (S/N) . By defining 

N 0 as the power spectral density where N 0 = N/BW, the detection 

index can also be written as d = 2t* (S/N q) . In Reference 2, the 

conditions required by the first model are called Case I and in 
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the following sections the first model is called the Case I 
model. A receiver that processes data such that it would 
implement a likelihood ratio decision rule under the conditions 
of the first model is called a crosscorrelation detector. If 
the description of the input noise is adequate, a Case I model 
can be used to obtain an estimate of an upper bound on a 
detection system's performance, since all the information 
necessary to define the signal is assumed to be known. 

Model II: In the second detection model, a sampled noise 

value is a value of an independent normal random variable with 
mean zero and with variance a*. And a sampled signal value is 
an independent random variable with mean zero but with variance 
(Jg. Thus, the input process corresponding to an observation 
consists of some number m of independent normal random 
variables Y 1 ,...,Y ro each with mean zero and each with variance 
a* when a signal is not present and with variance o' + Og when 
a signal is present. The result of applying a likelihood ratio 

decision rule in this model can be expressed in terms of a 

. . 2 

statistic X which is defined by X = S Y^. When a signal is 

not present, the statistic X/N has a chi-square distribution 

with m degrees of freedom. When a signal is present, the 

statistic X/(N+S) has a chi-square distribution with m 

degrees of freedom. So, in terms of these statistics, the 

conditional probabilities pf and p d are defined as follows: 

Pf- P(X^ > x*/N) and p d = P{X^ > (x*/N) [1/ (1+S/N) ] } 
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where is a chi-square random variable with m degrees of 

freedom, x* is a number which is determined by the decision rule 
and S/N is the signal-to-noise ratio. A receiver that would 
implement a likelihood ratio decision rule under the conditions 
of the second model is called an energy detector (or square law 
detector) . 

The mean of a chi-square random variable with m degrees of 
freedom is m and the variance is 2m. By the central limit 
theorem, as the number of degrees of freedom of a chi-square 
random variable becomes large it can be approximated by a normal 
random variable with the same mean and variance. Hence, for a 
sufficiently large sample size m or equivalently, for a 
sufficiently large time bandwidth product t-(BW), the 
conditional probabilities pf and p d can be approximated by 
p f = 1 - $(v*) and p d = 1 - ${ [1/ (1+S/N) ] (v* - d*)} where the 
threshold value v* = (x* - mo 1 ) / [ (2m) * a 1 ] and d = t* (BW) (S/N) 1 
if the concept of bandwidth is applicable. If, in addition, the 
noise N is significantly larger than the signal S, that is, if 
t* (BW) >> 1 and S/N << 1, then pf and p d can be 
approximated by: pf = 1 - *(v*) and p d = 1 - $ (v* - d^). 

In Reference 2, the conditions required for this approximation 
are called Case II and in the following sections the limiting 
form of the second model is called the Case II model. In 
Reference 2, an example is given that uses the Case II model to 
estimate the performance of a passive sonar system. 
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Model III: In the third detection model, a sampled noise 

value and a sampled signal value are values of independent random 
variables that are determined by independent Poisson processes 
observed for a time interval t. The noise process is 
characterized by a counting rate a, the signal process is 
characterized by a counting rate a s and the noise and signal 
processes are additive. Consequently, when the input is noise 
alone, the input is the value of a random variable Y that has a 
Poisson distribution with parameter at. And, when a target is 
not present, the input is a value of a random variable that has a 
Poisson distribution with parameter (a + a s )*t. For a 
likelihood ratio decision rule, pf = 1 - P(y*;at) and 
Pd = 1 - P[y*;(a + a s )*t] where P(y;0) symbolizes the 
cumulative Poisson distribution function with parameter 0 and 
y* is a threshold value of y determined by the decision rule. 

In terms of the above notation, both the mean and variance of a 
Poisson random variable are equal to 0. When 0 is large, the 
cumulative Poisson distribution function can be approximated by 
the cumulative distribution function of a normal random variable 
having the same mean and variance. Using this approximation for 
cases where at is large gives pj = 1 - #(v*) and 
Pd = 1 - #{[1/(1 + a s /a)*](v* - d^)> where v* = (y* - at)/(at) 5 
and d = a't 1 (a s /a) In addition, if the signal counting rate 
is significantly smaller than the noise counting rate, that is, 
if a*t » 1 and ag/a « 1, then pf and p^ can be 
approximated by: pf = 1 - #(v*) and pj = 1 - #(v* - d^). 
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This detection model might be used to describe a receiver whose 
input for an observation is the number of photons counted by a 
radiation detector in situations where at, the expected number 
of counts when no signal is present, is of the order of thirty or 
more. In this case, the target might be a wake segment for 
example . 

When a likelihood ratio decision rule is used in the three 
models discussed above, for the first model and under limiting 
conditions for the second and third models, the following result 
is obtained: pf = 1 - §(v*) and p<j = 1 - $(v* - d*) where the 

definition of v* depends on the noise power N for the first 
and second models. For a sonar receiver described by the first 
model, that is, by the Case I model: d = 2t* (BW) (S/N) . For a 

sonar receiver described under the limiting conditions for the 
second model, that is, by the Case II model, d = t • (BW) (S/N) 2 . 
So, in either a Case I model or a Case II model of a sonar 
receiver, the detection index d is a function of the time 
bandwidth product t* (BW) and signal-to-noise ratio S/N. Since 
sonar equations relate S/N to system, target and environmental 
parameters, a sonar equation can be used to relate S/N to these 
parameters in a model of a sonar receiver. 
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IV. General Detection Models 



The detection models that have been considered to this 
point are based on binary detection theory. After each 
observation, a receiver decides either that the input 
corresponding to the observation was noise or else it decides it 
was signal plus noise. However, in some detection systems this 
decision is delayed. In a computational sense, a model of such a 
detection system is generally complex relative to a binary 
detection model. To illustrate this, consider an active sonar 
system whose receiver includes an operator. Suppose the 
probability that the operator will detect a target echo has been 
determined in a laboratory experiment in which the operator was 
required to decide after each input corresponding to a resolution 
cell that either the input was a target echo (signal) and noise 
or the input was noise alone. In addition, suppose that under 
operational conditions the operator normally delays this 
decision. Then, in general, the probability that the operator 
will decide that the input corresponding to a resolution cell 
that contains a target is a target echo and noise will not be 
equal to the probability of the event in the forced choice 
experiment. And, in addition, the probability that the operator 
will decide the input corresponding to a resolution cell that 
does not contain a target is a target echo and noise will not be 
equal to the probability of this event in the forced choice 
experiment. Consequently, in general, the value of ! p u and 
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Pf for an operational environment will be different than that 
for the laboratory environment. 

One model that has been proposed to deal with this kind of 
situation defines the event that a receiver decides that the 
input corresponding to a resolution cell is signal and noise to 
be equivalent to the event that out of n consecutive 
observations at least k of them would result in the decision 
that the input was signal and noise in a forced choice 
experiment. The model is said to be based on an k-out-of-n 
detection criterion. With this criterion, the probability that a 
target will be first detected on the j th observation can be 
found as follows: Determine the 2^ sequences of forced choice 

responses that could result for a sequence of j consecutive 
observations. Next, determine the probability of occurrence for 
each sequence that first satisfies the k-out-of-n detection 
criterion on the j*-* 1 observation. The probability of first 
detection on the j th observation is equal to the sum of these 
probabilities. The cumulative probability of detection at the 
jth observation is the sum of the probabilities of first 
detection on the i^ observation for i = 1,2, ••• ,j. 
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V. Signal-to-Noise Ratio Detection Models 

In some radar and sonar detection models, for a 
specified value of p f , a minimum acceptable value of p^ is 
defined. This minimum acceptable value of p^ and the specified 
value of pf define what can be called a minimum acceptable 
signal-to-noise ratio (S/N) m if p^ is a nondecreasing 
function of S/N. In some sonar detection models, this 
signal-to-noise ratio in decibels is called the detection 
threshold DT. In symbols, DT = 10 log(S/N) m . If the minimum 
acceptable value of p^ is .5, then DT is usually called the 
recognition differential RD. The difference between the 
signal-to-noise ratio in decibels and RD (or DT) is called the 
signal excess SE. In symbols, SE = 10 log (S/N) - RD. 

One interpretation of signal excess is that for a 
localization region containing a target detection occurs with 
probability one if SE > 0 and with probability zero if SE < 0. 
This interpretation provides the basis for defining detection in 
the three encounter detection models that are discussed in the 
Section VII. A more consistent interpretation of signal excess 
is: If SE > 0, then the probability of detection p^ is 

greater than or equal to the minimum acceptable value (.5 if 
recognition differential RD is used to define signal excess) . 
For cases where p^ increases rapidly with signal excess in the 
neighborhood of zero signal excess, the two interpretations are 
essentially equivalent. For a discussion of this point as well 
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as a discussion of an operational case in which receiver 
decisions are delayed, see Reference 3. 

Signal excess ( signal-to-noise ratio) models provide a basis 
for general detection models, in particular, models that describe 
nonstationary noise and signal processes and randomly changing 
decision rules. This is illustrated by the models described in 
Section VII. In addition, signal excess models provide a basis 
for delayed receiver decision models. This is illustrated by the 
active sonar detection models in both Reference 4 and Reference 5 
that are based on a k-out-of-n detection criterion. In all of 
these models, the signal-to-noise ratio and the recognition 
differential are random variables. 

Using X(t) to represent a random variable corresponding to 
an index time t and a subscript to identify the random 
variable, for a passive sonar receiver, the signal-to-noise ratio 
in decibels associated with a decision at the index time can be 
expressed as follows: XsLtt) " x TL(t) " [X NL (t) - X DI (t)]. In 

this expression, SL represents source level, TL represents 
transmission loss, NL represents noise level and DI represen- 
ts directivity index. Since signal excess SE is defined to be 
the difference in decibels between the signal-to-noise ratio and 
the recognition differential (or detection threshold) , it too is 
a random variable and, for any decision time t, one can write: 

(1) X SE (t) = X SL (t) - X TL (t) - [X NL (t) - X DI ( t ) ] - X RD (t). 

The distributions of the random variables on the right side of 
Equation 1 determine the distribution of the signal excess. In 
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the passive sonar detection model described in Reference 6, 

XsL( t ) > x R[)(t) and, in effect, ^NL^) are normally 

distributed random variables while X TL (t) is a uniformly 
distributed random variable. In the three signal excess models 
that are described in Section VII, all of the random variables in 
Equation 1 are normally distributed. 

It is sometimes convenient to write Equation 1 as follows: 

(2) X SE (t) = SE ( t ) + X ( t ) 

In Equation 2, SE(t) is the expected value of the signal excess 
determined by the following expected value equation: 

(3) SE ( t ) = SL(t) - TL(t) - [NL(t) - DI (t) ] - RD(t) 

where each term on the right represents the expected value of the 

indicated random variable and X(t) is a random variable that 

determines the stochastic character of the signal excess. Since 

SE(t) is the mean of Xg E (t); by Equation 2, the mean of X(t) 

is equal to zero and the standard deviation of X(t) is equal to 

the standard deviation of Xgg(t). If a represents the 

standard deviation of Xg E (t) and the random variables on the 

right side of Equation 1 are statistically independent, then 
2 2 2 2 2 2 

a = °SL + a TL + a NL + a DI + a RD- This relation has been used to 
determine a standard deviation for the signal excess in 
operational models. 
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VI. General Encounter Models 



A basic problem associated with search modeling is that 
of determining the probability that a target will be detected by 
at least one detection system in an encounter with one or more 
detection systems on or before the n^h decision in the encounter 
has occurred regarding the presence or absence of a target in a 
localization region. In the encounter models that are considered 
in this report, during a search, observations are made of a 
series of localization regions. For localization regions that do 
not contain a target, the signal observed during the observation 
of the region is zero. In addition, in these models, the time to 
resolve a false alarm is zero. However, p d and pf are 
considered to be determined by some criterion such that Pf is 
less than one. Consequently, although the time to resolve a 
false alarm is zero, the cost associated with a false alarm is 
not zero. 

Using an order number rather than a time to index a decision 
and a random variable N to represent the decision order number 
at which detection first occurs, the cumulative probability of 
detection for an encounter is given by: P(N < n) = E P(N = i) 
where the sum index i = 1,2, ••• ,n. It also given by: 

P(N < n) = 1 - (1 - gi )(l - g 2 ) ... (1 - g n ) where 
gi = P(N = i|N < i-1) is the probability of the event detection 
on the i*-* 1 decision conditioned on the event no detection on an 
earlier observation. The second expression for P(N < n) is 
generally of greater interest than the first expression, since 
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can usually be more directly related to operational 
parameters such as target range and target environment that 
determine a targets detectability than can P(N = i) . If a 
target is not present in the localization region that corresponds 
to the i^h decision, then g^ will be zero. The second 
expression above for P(N < n) can also be written in the 
following form: P(N < n) = 1 - exp[E ln(l - gjJ ] where the sum 

index i = 1,2, ••• ,m. In this expression, ln(l - g^) can be 

approximated by -g^ for g^ << 1. With this approximation, 

P(N < n) = 1 - exp(-E g^) . This expression provides a basis for 
developing a continuous analog for P(N < n) that has been used 
to describe detection systems for which decisions can be 
considered to occur continuously. That is, for cases where the 
time of the observations corresponding to a decision and the time 
between decisions are both negligible relative to the average 
time of an encounter. The analog can be developed as follows: 
First, replace the decision index by a time index. Next, let St 
be the time between decisions. Then, the cumulative probability 
of detection after the n th decision can be expressed by: 

P(N < n) = P (T < t) where T = N*6t and t = n*5t. Next, 
define a detection rate function (a probability of detection per 
unit time) by the expression: (l/St)P[T = t|T < (n-l)*St). 

If T is considered to be a continuous random variable, then 
this expression can be used as a guide for defining a continuous 
detection rate t ( t) . The definition is: 
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(4) 



T ( t ) = lim{ (l/6t)P(t < T < t+6 t | T < t) } 



where the limit is for 6t approaching zero and P(T < t) is 

such that P(T < t) = P(N < n) when t corresponds to n for 
all values of n that are involved and P(T < t) is a 
nondecreasing function of t. Equation 4 implies the following 
differential equation: dp(t)/dt = [1 - p(t)]T(t) where 

p(t) = P (T < t) . The solution to this differential equation is: 



where t is the time index for a decision during an encounter, to 
and tl are times of the encounter and tl > to. A T(t) that is 
based on a visual detection model is described in Reference 7. 

If the detection capability of a detection system is assumed to 
depend on a target's position relative to the detection system 
during an encounter but not to depend on the clock time, then the 
time index of a decision can be a relative index that determines 
a target position that is associated with a decision rather than 
the clock time associated with the decision. 

The above results apply to the case of an encounter between 
a target and a detection system consisting of more than one 
sensor. However, if the sensors are not located on a single 
platform or at a single point, it is generally more convenient to 
describe encounters of this kind in terms of the encounters 
between the target and the individual sensors. In either case, 
if the event that a sensor detects a target during an encounter 
is not independent of the this event for other sensors, in order 



(5) 



P (T < tl) = 1 - P (T < tO) 
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to describe this in an encounter model, the correlation between 
this sensor and the other sensors must be specified. This has 
been done in some models as follows: First determine the 

probability of detection during the encounter for each sensor 
involved. Then determine the probability that at least one of 
the sensors will detect the target during the encounter for the 
case where the random factors that determine detection are 
independent and again for the case where they are completely 
dependent. Let be the probability that the it* 1 sensor 

detects the target during the encounter and let the probability 
that at least one sensor detects the target be Pj in the first 
case and P^ in the second case. For the first case, 

Pj = 1 - (1 - Pj) (1 - P 2 ) ••• (1- P n ) where n is the number 

of sensors involved. And, for the second case, 

P D - P m where P m > P^ for i = 1,2, ••• ,n. Next, combine the 

two probabilities as follows: P = a*P D + (1 - a)*Pj. A way to 

determine a value for a is described in Reference 8. 
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VII. Three Signal Excess Encounter Models 

In the three models described in this section, 
detection is defined in terms of signal excess as it is in 
Section V. Each of the models determines a cumulative 
probability of detection for a target in an encounter with a 
passive sonar system. An observation in the models is indexed by 
time and the index can usually be considered to be the time at 
the end of the observation. During an encounter, observations 
are made of a series of localization regions. For localization 
regions that do not contain a target, the signal observed during 
the observation of the region is zero. The time to resolve a 
false alarm is zero. However, false alarms are not ignored in 
that the value of RD (or DT) is determined by some specified 
false alarm probability. Consequently, although the time to 
resolve a false alarm is zero, the cost associated with a false 
alarm is not zero. In determining signal excess in the models, 
it is convenient to use Equation 2. For each decision in an 
encounter, there is a random variable X(t) defined by Equation 
2 that determines the random character of the signal excess. 

For a sequence of decisions, the set of these random variables 
ordered by their time index constitutes a stochastic process. 

The joint distributions of these random variables determines the 
nature of this process. In the three encounter detection models 
described in this section, the process is called a lambda sigma 
jump process. The jumps in the series occur at times determined 
by a Poisson process with a mean rate lambda. This implies that 
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the time between jumps is a random variable with an exponential 
distribution and that the expected times between jumps t is 
equal to the reciprocal of lambda. The time series that are 
generated by lambda sigma jump processes are represented by 
Figure 2 below. 




Figure 2. A time series representing a realization of a lambda 
sigma jump process. On the plot, o in dB equals one unit on 
the vertical axis and i equals one time unit on the horizontal 
axis . 



From Figure 2, note that the observed values of neighboring 
random variables are equal unless a jump has occurred between 
them. When a jump occurs, the first random variable after the 
jump is normally distributed with mean zero and variance a 2 and 
it is independent of all the random variables before the jump. 
Conditioned on a jump pattern, this random variable and a] 1 the 
random variables between it and the next jump are dependent and 
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the correlation coefficient between any pair is one. That is, if 
the value of the signal excess is known at some time, then all of 
the values between the last jump before that time and the first 
jump after that time are also known. However, since the jumps 
occur randomly, knowing the value of the signal excess with 
certainty at some time does not determine the values of the 
signal excess with certainty at neighboring times. In the 
unconditioned case, the correlation coefficient between the 
random variables X(t) and X(t+T) is equal to 1/e. For this 
reason, t is referred to as a relaxation time. If only the 
transmission loss is considered to be determined by a stochastic 
process, then only the transmission loss will be a random 
quantity and ^ is the expected time between jumps in the random 
component of the transmission loss. It appears that the use of 
the lambda sigma jump process is based more on past practice than 
on experimental justification. In this regard, see Reference 9. 
By referring to Equation 1, it can be seen that the lambda sigma 
jump stochastic process is equal to the sum of the stochastic 
processes that determine the random variables on the right side 
of this equation. But, the nature of the processes should depend 
on the encounter that is being modeled. This suggests that an 
adequate description of an encounter using a lambda sigma jump 
process might not be possible in all cases. 

In the three encounter models described below, detection is 
defined in terms of signal excess as described in Section IV and 
decisions are indexed by a time that can usually be considered to 
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be the time of the decision. During an encounter, observations 
are made of a series of localization regions. For localization 
regions that do not contain a target, the signal observed during 
an observation of the region is zero. In the models, the time to 
resolve a false alarm is zero. However, false alarms are not 
ignored in that the value of RD (or DT) is determined by some 
specified false alarm probability. Consequently, although the 
time to resolve a false alarm is zero, the cost associated with a 
false alarm is not zero. 

The First Passive Sonar Encounter Detection Model: This 

model describes an encounter as a series of decisions. In the 
model, the integration time that determines the recognition 
differential is equal to the duration of the observations 
corresponding to a decision. For an encounter in which SE(t) 
is unimodal and in which the time of the single maximum is prior 
to the end of the encounter, it is shown in Reference 10 that 
the probability that detection will occur during the encounter 
is given by the following equation: 

(6) p = 1 - [ (1 - p c )/(l - p- Pc )](l - P • Pl ) ••• (1 - P-p m ) 
where P = 1 - exp(-St/T) and p^ = $[SE(ti)/o] for 1 = 1,2,3, 
,m. Here, St indicates the duration of an observation and 
$ indicates the standard normal cumulative distribution function 
as before. The integer c is the index of the decision time t c 
for which SE(t c ) is greater than or equal to SEft^) for any 
time t^ and t c < t m . 
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As t approaches zero, 3 approaches one and Equation 6 
approaches this form: 

(7) p = 1 - (1 - Pi ) ••• (1 - p m ) . 

In this limit, the signal excess random variables are all 
independent. Note that Equation 7 applies without the condition 
that SE ( t) be unimodal. 

As t approaches infinity, 3 approaches zero and 
Equation 6 approaches this form: 

(8) P = P c - 

In this limit, the correlation coefficient between any pair of 
signal excess random variables is equal to one. Note that 
Equation 8 applies without the condition that SE(t) be 
unimodal. Equation 8 is said to define a complete dependence 
encounter model. 

The Second Passive Sonar Encounter Detection Model: This 

model is in a sense a third limiting form of the first passive 
sonar encounter detection model. In this limit, the time between 
decisions approaches zero. However, in this limit the 
integration time that determines the recognition differential is 
not equal to St and it does not approach zero. It is, in 
effect, chosen by the user of the model through the user's choice 
of the value for the recognition differential. 

For an encounter that begins at tl and ends at t2 and 
for which SE(t) is unimodal, it is shown in Reference 10 that 
for this limit. Equation 4 has the following form: 



(9) p = 1 - [1 - p(t c )] exp[- ( 1 / t ) _ p(t) dt] 
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where p(t) = $[SE(t)/o] and where now t c is the encounter 
time such that SE(t c ) is greater than equal to SE(t) for any 
other encounter time t and t2 > t c . 

The Third Passive Sonar Encounter Detection Model: This 

model describes an encounter between a target and a passive sonar 
detection system in which detection occurs at a time t if the 
average value of the square of the continuously observed 
signal-to-noise ratio over a time interval of length u from 
time t-u to time t is greater than or equal to the square of 
the signal-to-noise ratio that determines the recognition 
differential for that time interval. With R(s) the random 
signal-to-noise ratio at a time s and R m (u) the random 
signal-to-noise ratio that determines the random recognition 
differential for the time interval u (the integration time) , 
detection at time t is the following event: 

(10) (1/u) J t _ u [R(s)/R(u)m] 1 ds > 1. 

The random integrand in the expression that defines the event is 
related to the random signal excess at the time s for an 
integration time u. The relation is: 

(11) 10 log [R(s)/R m (u) ] * = 2[SE(s;u) + X(s)] 

where SE(s;u) is the expected value of the signal excess at a 
time s for an integration time u and X(s) is the random 

component of the signal excess at the time s. In the model, 

X(s) is determined by a lambda sigma jump process and SE(s;u) 
is determined by an expected value sonar equation with a 
recognition differential: RD(u) = 10 log r m (u). Here, r m (u) 
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is the expected value of the signal-to-noise ratio that gives a 
probability of detection equal to .5 for an integration time u 
and a specified probability of false alarm pf. With the signal 
detection process described by a Case II signal detection model, 
the detection index necessary to give the required operating 
point (pf,.5) is related to the integration time u and the 
signal-to-noise ratio r m (u) by: 

(12) d = u- (BW) [ r m (u) ] 1 

where BW is the bandwidth of the receiver. For a spectrum 
analyzer, BW would be the bandwidth corresponding to a given 
frequency resolution and d would be the detection index 
required in order to be at the operating point (pf,.5) for a 
signal that was contained within a bandwidth BW. Since d is 
fixed, because of Equation 12: 

(13) RD (u) = 5 log (t 0 /u) + RD(t 0 ) 

where t 0 is a reference integration time. By using Equation 13, 
Equation 10 becomes: 

^ (1/5) [X(s) + SE(s;t 0 ) - 5 log (t 0 ) ] 

(14) J t _ u 10 ds > 1 

after moving the time origin so that t > 0. In this equation, 
u = t for t < t 0 and u = t 0 for t > t 0 and SE(s;t 0 ) is 
the expected value of the signal excess at the time s for an 
integration time tg. In an encounter, detection occurs the 
first time that Relation 14 is satisfied. 

As is pointed out in Reference 11, the appeal of the Third 
Passive Sonar Encounter Model relative to the Second and First 
Passive Sonar Encounter Detection Models is that it appears to 
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more closely describe the detection process in passive sonar 
detection systems that display their processed data to an 
operator in a continuous manner. However, results reported in 
Reference 12 indicate that the difference between the three 
models may not be significant in some types of encounters. 
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VIII. Straight Line Encounters 

Suppose a target's detectability depends on its range 
from a detection system and that the probability of detection is 
effectively zero beyond a range r m for any target azimuth. In 
this report, an encounter between the target and the detection 
system is the event that the range between the target and the 
detection system is less than or equal to r m . In addition, 
suppose r m is small enough so that when the target and the 
detection system are having an encounter they can be considered 
to be moving on planes parallel to a tangent plane to the earth's 
surface at some point in their vicinity. If this is the case, 
then while the target and detection system maintain a constant 
course and speed during an encounter, the encounter is called a 
straight line encounter. 

A straight line encounter can be described in terms of a two 
dimensional rectangular coordinate system whose plane is parallel 
to the tangent plane to the earth. If the coordinate system is 
stationary relative to the detection system with the detection 
system located at the origin and is oriented so that the target's 
motion is parallel to the y-axis and is in the positive 
y-direction, then the target's x-coordinate during a straight 
line encounter will be constant. The constant is equal to the 
target's horizontal range at the closest point of approach (CPA) 
on the straight line track on which the target is moving relative 
to the detection system during the encounter. This range is 
called the target's lateral range. 
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A complete straight line encounter is one in which a target 



begins an encounter on a straight track at a point on the 
relative track whose range from a detection system is equal to 
r m , continues to CPA and then continues to a point on the 
relative track whose range from the detection system is again 
equal to r m . Let p(x) be the cumulative probability that a 
target is detected by a detection system in a complete straight 
line encounter in which the target's lateral range is x. Then 
the function p(x) defines what is called a lateral range curve 
or lateral range function. 

Let p be the probability that a target is detected during a 
complete straight line encounter. If the lateral range of a 
target in a straight line encounter is assumed to be a continuous 
random variable X with a uniform distribution that extends at 
least from -r m to r m where r m is the maximum detection 
range, then the probability that a target will be detected during 
a complete straight line encounter is: 



where the limits of integration can be used since the value of 
p(x) is zero for x < -r m and for x > r m . The quantity: 



is called the sweep width. Using the definition of W given by 
Equation 16, the probability that a target will be detected 
during a complete straight line encounter under the conditions 
that determine Equation 15 is p = (l/2r m ) *W. 



(15) 




(16) 
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IX. Two Intermittent Signal Encounter Models 

In the intermittent signal encounter models that are 
described in this section, an encounter is a complete straight 
line encounter and during an encounter a target emits a signal at 
various times. Two cases are considered: In the first case, the 
signals are emitted periodically, the signals are of length St 
and the time between signals is t where t > St. In the second 
case, the signals are instantaneous and the signals are emitted 
at times determined by a Poisson process for which the expected 
time between signals is equal to t. In the model, the 
detectability of a target signal depends on a target's horizontal 
range from a detection system, but on no other factors. If a 
signal is emitted while the target is within a range r, it will 
be detected. For a continuous signal, the lateral range function 
of a detection system for a target is: P(x) = 1 for |x| < r 
and p(x) =0 for |x| > r where the horizontal range r is 
determined by the characteristics of the detection system and the 
target. The geometry for an encounter is shown in Figure 1 
below. 

For an intermittent signal, the length of a target's track 
relative to a detection system where if a signal is emitted it 
will be detected is 2(r l - x 1 )* + w*St with w the speed of 
the target relative to the detection system. So, a target's 
exposure time during an encounter is (2/w) (r 1 - x 1 )^ + St. 

For a periodic signal, there are two cases. In the first 
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case, r > w(t - 6t)/2. In this case, the exposure results in 
the following lateral range function: 
p(x) = 0 for |x| > r 

(17) p ( x ) = 1 for | x | < (r 1 - [ w • ( t - &t)/2]*)i 

p(x) = (2/w*T)(r J - x 1 )’ + St/T otherwise 
In the second case, r < w*(t - 5t)/2. In this case, the middle 
equality in Equation 17 does not apply. 

For a signal that is instantaneous and whose emission times 
are determined by a Poisson process, the exposure results in the 
following lateral range function: 

p ( x) = 1 - exp{ - [ ( 2/w t ) (r 1 - x 1 )*]} for |x| < r 



( 18 ) 

p(x) = 0 for | x | > r 

for all values of r. A case in which the signal length St is 
greater than zero is discussed below. 




Figure 3 . The encounter geometry for the two intermittent signal 
models described above. 
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For signals whose emission times are determined by a Poisson 
process and whose length is fit, signals can overlap. If this is 
allowed, then Equation 18 can be modified to describe this case 
by adding fit/T to the term in the exponent of Equation 18 that 
is within the square brackets. In particular, note that this 
modified Equation 18 can be approximated by the bottom equality 
in Equation 17 when (2/w*T)(r l - x’)* + St/T << 1. This 
implies that when the expected time t between signals is large 
relative to the exposure time (2/w) (r l - x 1 )* + fit the periodic 
emission model and the random emission model are effectively 
equivalent. 
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X. A Random Search Model 



A search of a region in which a target moves on a track 
that consists of a number of straight segments placed in such a 
way that in a limiting sense every section of the region is 
equally likely to be searched on a segment is referred to as a 
random search in Reference 6. Representation of a search region 
with the track segments that could be imagined to be the tracks 
of a random search are shown in Figure 4 . 




Figure 4. A search region and a track that could be described as 
a random search track. 

Two developments of a model to describe this kind of search 
are contained in this section. The first development is based on 
the following conditions: 1. A target is at a fixed position 

within a defined search region. 2. A searcher's track is a 
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sequence of straight line segments that are within the search 
region. 3. The searcher's detection system is such that while on 
a track segment, a rectangle is searched that is contained within 
the search region, is of length equal to the length of the track 
segment and is oriented so that its long axis is parallel to the 
track segment. 4. The probability that the searcher's detection 
system will detect a target while on a track segment with a 
search rectangle that does not include the target is zero. The 
probability that the searcher's detection system will detect a 
target while on a track segment with a search rectangle that 
includes the target is p(x) where x is the target's lateral 
range for the track segment and p(x) is the lateral range curve 
for a complete straight line encounter with this lateral range. A 
representation of a search rectangle is shown in Figure 5 below. 
5. The track segments are located in such a way that the event 
that the target is within the search rectangle associated with a 
track segment is independent of the event that the target in the 
search rectangle associated with any other track segment. And 
the probability of the event is equal to the ratio of the area of 
the search rectangle to the area of the search region and, given 
a target is within a search rectangle, its position is uniformly 
distributed over the rectangle. 

Condition 4 implies that the random search model is based 
on the concept of a complete straight line encounter. The 
definition of an encounter that is intended here is that given in 
Section VI. This implies that in the random search model the 
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time to resolve a false alarm is zero. However, for the model, 
p <3 and pf are considered to be determined by some criterion 
such that pf is less than one. Consequently, although the time 
to resolve a false alarm is zero in the model, the cost 
associated with a false alarm is not zero. (A concise model that 
accounts for the time to resolve false alarms is described in 
Reference 13.) Condition 4 also implies that when a searcher is 
on a track segment with a search rectangle that contains a 
target, the encounter is a complete straight line encounter. And 
Condition 5, which can be considered to specify a random 
arrangement of the track segments, implies that when this is the 
case, for the complete straight line encounter, the target's 
lateral range is a random variable that is uniformly distributed 
between -b/2 and b/2 where b is width of the search 
rectangle (the dimension of the rectangle perpendicular to the 
associated track segment) . 



Swept Area 



Track 

Figure 5. A track segment and its associated search rectangle, 
that could correspond to a search with an aircraft mounted 
infrared detection system. 
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Based on the above considerations, the probability that a 
target will be detected while a searcher is on a track segment 
with an associated search rectangle that contains the target is 
given by: 



pOO 



(19) J — oo P (x)f x (x)dx = W/b 

where fx( x ) = l/ 2 for -b/2 < x < b/2 and f x (x) = 0 and 
p(x) = 0 otherwise. Note that the left side of Equation 19 
applies to any complete straight line encounter in which the 
target's lateral range for the encounter is considered to be a 
random variable with a distribution determined by the probability 
density function f x (x) • If it is not given that the target is 
within the search rectangle associated with a track segment, then 
the unconditional probability that the target will be detected on 
the track segment is given by: (W/b) (SA/A) where & A is the 

area of the search rectangle associated with the track segment 
and A is the area of the search region. With 1 the length of 
the rectangle, 6A = b*l and the probability becomes: (W*l)/A. 

Then, since the event that the target will be in the search 
rectangle of a track segment is independent of the event that it 
will be in the search rectangle of any other track segment, the 
probability p that a random search consisting of m track 
segments will detect the target is given by: 



1 - [1 - (W-l^/AHl - (W*l 2 )/A] ••• [1 - (W-l n )/A] where 1* 
is the length of the i th track segment. The probability is 
also given by: p = 1 - exp{E ln[l - (W*l^)/A]). If for 
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i = 1,2, ••• , n (W-li)/A << 1, then this expression can be 

approximated by: 

(20) p = 1 - exp[-(W-l)/A] 

where the sum index i = 1,2, ••• ,m and 1 = Z lj_ is the track 
length of the search. Equation 20 is known as the random search 
formula. 

The second development is based on the following condition: 
The detection rate of a searcher's detection system for a target 
is given by: T(t) = [W*v(t)]/A where t (t) is the continuous 

detection rate function in Equation 5. Using the value for t (t) 
and Equation 5, the random search formula given by Equation 18 is 
easily found. 

The advantage of the method used in the second development 
of the random search formula over the method used in the first 
development is that, in general, it can be more easily used to 
determine the probability of detection for a search that can be 
defined in terms of a detection rate. Reference 14 contains an 
example of an application of Equation 5 to a case where the 
search region expands with time. 
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XI. Ladder and Barrier Search Models 



In some barrier searches, the barrier search track is a 
ladder search track relative to a reference system that moves 
with the target. This fact is used in the barrier search model 
development that follows the two ladder search model developments 
below. The first ladder search model is referred to as an ideal 
ladder search model because of the idealizations that are 
involved in its description of a ladder search. The second 
ladder search model is referred to as a degraded ladder search. 

It can be considered to describe a ladder search track in which 
navigational errors result in omissions and overlaps in coverage. 

An Ideal Ladder Search Model: The model is based on the 

following conditions: 1. A ladder search region is a rectangle 

that contains a fixed target. 2. During a search of the region, 
the searcher's detection system moves on a set of m parallel 
track segments of length b separated by a distance s. 3 . As 
the detection system moves along a track segment, it searches a 
rectangular strip of length b and width s within the search 
region. 4. The m rectangular strips that correspond to the m 
track segments completely cover the ladder search region with no 
overlap. 5. If a target is within the rectangular strip 
corresponding to a track segment, then there will be a complete 
straight line encounter between the target and the detection 
system when the detection system moves along the track segment 
and the lateral range of the encounter will be uniformly 
distributed across the width of the strip. If the target is not 
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in the rectangular strip, then there will not be an encounter and 
the probability that the target will be detected while the 
detection system is on the track segment is zero. 




— s — 

Figure 6. A schematic representation of a ladder search geometry 
for a case in which the ladder search track segments are 
superimposed on and bisect their corresponding rectangular 
strips . 

Since targets outside of the rectangular strip that 
corresponds to a track segment cannot be detected while a 
detection system is on the track segment because of Condition 5, 
in the model, the sweep width W of a searcher's detection 
system must satisfy the relation W < s . In particular, W = s 
only holds when the detection system detects a target that is in 
a rectangular strip with probability one for any target lateral 
range. This kind of detection system is sometimes referred to as 
a cookie cutter detection system. However, this terminology can 
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be misleading since it suggests the detection system detects 
equally well for all azimuths. But this is not a requirement on 
the system in order that W = s. 

The ideal ladder search model implies that if the conditions 
of the model are satisfied, then the probability p that a 
target will be detected by a an ideal ladder search is given by: 

(21) p = W/s 

where W/s < 1. The quantity W/s is called the coverage 
factor in this case. 

A Degraded Ladder Search Model: The above model implies 

perfect navigation in addition to other idealizations. A model 
of a ladder search is given in Reference 6 that could be used 
for cases in which this is a poor assumption. The model which is 
referred to here as a degraded ladder search model can be 
considered to describe navigational inaccuracies in terms of 
omissions and overlaps of the rectangular strips. It can be 
developed as follows: Consider a random search in the ladder 

search region whose track length is equal to the search track 
length required to complete an ideal ladder search, that is, a 
track length 1 = m*b. The degraded ladder search model 
describes the result of omissions and overlaps in a ladder search 
to be such that the probability of detection for this random 
search is equal to the probability of detection for the degraded 
ladder search. Consequently, since the area of the ladder search 
region is m*s*b, for the degraded ladder search model: 

(22) p = 1 - exp(-W/s). 
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Here, the requirement that the coverage factor W/s < 1 for 
Equation 21 can be relaxed. However, it should still be 
considered as an approximate condition. 

The condition that the target be fixed within the 
rectangular search region is critical to both Equation 21 and 
Equation 22. However, these results are also applicable to a 
search for a moving target under the conditions that are 
described next. 

A Barrier Search Model: A target moves with a constant 

course and a constant speed u. Both the target's course and the 
target's speed are known by a searcher. The searcher establishes 
a barrier of width b that is perpendicular to the target's 
track and moves on the barrier with a speed v > u. The barrier 
is designed so that in a reference system relative to the target 
the barrier search is a ladder search that satisfies the 
conditions for a ladder search that are given above. There are 
two cases to consider: 1. The barrier is established in front of 
the target. 2. The barrier is established behind the target. 

From the search geometry for a barrier established in front 
of the target, it can be seen from Figure 7 below that 
9 = sin -1 (u/v) and d = v*t where t = s/ (v + u) is the time to 

move from one search leg to the next. The angle 0 and the 

perpendicular distance d that depend on u, v and s, and the 

width of the barrier b are the quantities that are required in 
order to establish the barrier operationally. 
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b - 

Figure 7. A barrier search track shown for a barrier established 
in front of the target. The track is shown in a reference system 
fixed relative to the earth. 

For a barrier that is established in front of a target, one 
of three barrier types will result. A barrier's type is 
determined by the relation of the distance d to the distance 

g = ut where the time t = b/(v s - u 1 )* is the time to complete 

a search leg (cross the barrier) . The barrier type is determined 
as follows: 1. For g < d, the barrier is an advancing barrier. 

2. For g = d, the barrier is a stationary barrier. 

3. For g > d, the barrier is a retreating barrier. 

For a barrier established behind the target, there is only 
one barrier type and it is called an overtaking barrier. For an 
overtaking barrier, 8 = sin -1 (u/v) as for a barrier established 
in front of the target. But, for an overtaking barrier, 
t = s/(v - u) and d = vs/(v - u) . 
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Given the target crosses the barrier, the probability of 
detection for an ideal barrier search is given by Equation 21 
and the probability for a degraded barrier search is given by 
Equation 22 where the terminology refers to the nature of the 
ladder search in the reference system moving with the target. 

A discussion of an application of these two equations to a search 
for a magnetic anomaly target is given in Reference 15. 
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XII. A Target State Estimation Procedure 

A target state estimation procedure based on bearing 
observations is developed in this section that generates point 
estimates of a target's position and velocity vector coordinates 
in a rectangular coordinate system. The procedure is based on a 
model in which bearing errors are unknown and are not determined 
by random variables with known distributions. Because of this, 
confidence regions for the estimates are not generated by the 
procedure. However, for a moving target, it illustrates some 
general characteristics of bearings only target motion analysis 
(TMA) . 

The target state estimation procedure model is defined by: 

1. The target moves in a plane with a constant but unknown course 
and speed. 2. Observations of the target are made from known 
positions at known times. 3. The observations provide only 
target bearings with unknown errors. 4. Target position 
estimates [x t (i) ,y t (i) ] for i = 1,2, ••• , n and target 

velocity component estimates u x and Uy are chosen so that the 
sum S = 2 d{ is a minimum where d^ is the algebraic distance 
between the estimated position and the bearing line corresponding 
to the observed bearing. 

Because of the requirement that the target move with 
constant course and speed during the encounter, the number of 
independent estimates is reduced from 2m to 4: x-t-(l), y^fl), 

u x and Uy. These four estimates constitute a set of target 
state parameters that define a target's motion by the relations: 
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x t(i) = xtfi) + u x*( fc i " t l > and yt(i) = Yt( 1 ) + u y* ( fc i " fc i) • 
The geometry on which the procedure is based is shown in Figure 8 
where d A = [x t (i) - x Q (i)]-cos ©i - [y t (i) - y 0 (i)]-sin 0^. 

North 




Figure 8. The geometry of the target motion analysis model. 

To determine estimates of the target state parameters, take 
the partial derivative of S with respect to each of them. Then 
set the four partial derivatives equal to zero. This creates 
four linear equations in the four unknown estimates. In matrix 
notation, the equations are: AX = B where the elements of X 

are x t (l), y t (l), u x> and u y A necessary condition for a 
unique solution for X is that m > 4. Otherwise, the 
determinant of A will be equal to zero. 

Now, suppose the observations are at positions and times 
that correspond to the positions and times of an observer moving 
on some constant course at some constant speed (including zero 
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speed). In this case, the observation position coordinates are 
related by the following equations: x Q (i) = x Q (l) + v x (t^ - t^) 

and y Q (i) = y Q (l) + v y( fc i “ fc i) where v x and v y are the 
required velocity components of the observer. Using these 
equations of motion, the matrix equation AX = B can be 
transformed to the matrix equation Ax'= 0 where the elements 
X £2 of the matrix X are related to the elements of the matrix 
A* by the equations: x^ = x t (l) - x Q (l), X 21 = yt(l) “ y 0 ( 1 )/ 

X 31 = u x " v x and X 41 “ u y ” v y 

Since the equations represented by Ax' = 0 are homogenous, 

the equations do not have unique solutions and consequently the 
equations represented by AX = B do not either. However, if 
there is at least one observation whose time and position is not 
determined by the above equations of motion, then the 
transformation from X to x' cannot be made, and in general a 
unique solution for X can be found. If the observations are 
made from a platform that is moving with a constant course and 
speed, this condition can be achieved by either changing the 
course, the speed or both prior to completing the observations. 

Estimation models that describe bearing error in terms of a 
random variable with a specified distribution provide a basis for 
determining a confidence region to associate with the point 
estimates. A model is developed in Reference 16 that does this, 
however, it requires that either observations of a target bearing 
be made from two or more locations simultaneously or that the 
target be stationary relative to the observation points. 
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XIII. Position Distributions That Change with Motion 

In this section, target motion models are discussed 
that provide a basis for determining position distributions that 
change with time. In the first model, a target moves in a plane 
with a constant course and speed. In the second model, a target 
moves in a plane with a constant course and speed during a time 
step and at the end of the time step changes either its course or 
speed or both. 

Position Distributions based on the First Motion Model : Two 

cases of this model are considered that provide the basis for 
determining a position distribution analytically. In both cases, 
the target's position coordinates in a rectangular coordinate 
system are determined by a circular normal distribution at time 
zero. And, in both cases, during the motion the target's 
velocity does not change. 

The First Case: In this case, the targets 's coordinates 

X(t) and Y(t) at time t are determined by the vector sum 
with components: X(t) = X(0) + U x *t and Y(t) = Y(0) + Uy*t 

where X(0) and Y(0) are the target's coordinates at time zero 
and U x and Uy are the target's velocity components. The 
coordinates X(0) and Y(0) are determined by a circular normal 
distribution whose mean vector is (0,0) and whose standard 
deviation is a. And the velocity components U x and Uy are 
determined by a circular normal distribution whose mean vector is 
(u x ,Uy) and whose standard deviation is o u . The random vectors 
[X(0),Y(0)] and (U x *t,Uy*t) are independent and both are 
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normal random vectors. Since they are normal random vectors, the 
random vector [X(t),Y(t)] is also a normal random vector. And, 
since they are independent, its characteristic function is equal 
to the product of their characteristic functions. This implies 
that [X(t),Y(t)] has a circular normal distribution whose mean 
vector is (u x *t,Uy*t) and whose variance is (o’ + a^) . For a 
more complete discussion of the basis for the above argument, see 
Reference 7 . 

The Second Case: The procedure that is used to determine 

the distribution of the random vector [X(t),Y(t)] for the 
second case can also be used to determine it for the first case. 
It is based on the following consideration: In general, the 

joint density function of a target's coordinates at some time t 
greater than a zero of time is determined as follows: 



poo poo 



( 23 ) f X(t) , (t) (x,y;t) = J-ooJ-oofx(O) ,Y(0) ( r /s;°) f Ux ,uy ( v ' w ) dvdw 

where r = x - v*t and s = y - w*t with v = u x and w = Uy. 
Equation 23 can be interpreted as follows: After multiplying 

the right-hand side of Equation 23 by SrSs, in the limit, it is 
the sum over all pairs of values (v,w) of the product of the 
probability that the target's velocity components (v,w) are in 
an element SvSw and the probability that the target's position 
coordinates (r,s) at time zero are in an element of area SrSs. 
This is the position such that the targets 's position coordinates 
(x,y) at time t will be in an element of area 6xSy given 
that the target's velocity is (v,w). And, in the limit, this is 
equal to the left side of Equation 23 after multiplying the left 
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side by 6x6y. But, since 6r6s is equal to 6x6y in the 
limit, Equation 23 results. In effect, 6r5s is translated 
from (r,s) to (x,y) without rotation or distortion. 

In the first case, the joint distribution of the of the 
velocity components (U x ,Uy) is circular normal with mean vector 
(u x ,Uy) and standard deviation o u and the joint distribution 
of [ X ( 0 ) , Y ( 0 ) ] is circular normal with mean vector (0,0) and 
standard deviation o. In the second case, the target's speed 

u is known but its course is equally likely to have any value 

between 0 and 2 tt, so U x = u*sin $ and Uy = u*cos $ where * 

is a uniform random variable and 0 < i < 2ir . In this case, only 

the single random variable $ is required to determine the 
target's velocity. The integral of Equation 23 in polar 
coordinates is a single integral over $ and the integrand is 
(1/2TTO*) * exp [ — ( r * + s * )/2o* ] • (l/2ir) where r = x - u*t*sin <p 
and s = y - u*t*cos <p. The result of the integration is: 

(24) ( l/2iro 1 ) • exp{ - [x* + y* + (u • t) 1 ]/2o J ) I 0 [ (x * + y’J^u^t/a*] 

where I 0 indicates the hyperbolic Bessel function of zeroth 
order. Expression 24 is plotted in Reference 6 for several 
values of t. There, a single variable is substituted for the 
distance of the target from the origin (x* + y*)^ . It is 
measured in units of o and the time is measured in units of 
o/u. The plot of the expression shows a limiting characteristic 
of the distribution that can be indicated as follows: Multiply 

and divide the resulting expression by exp(-z*U‘t/a’) where 
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z = ( x a + y a )^, the target's range from the origin. The result 
of the operation can be expressed in the following form: 

(25) 1/ ( 27 ra a ) exp{ - [ 1/ (2a 1 ) ] ( z - u • t) * ) I 0 ( z • u • t/o a ) exp (-r *u • t/o a ) 

where the second factor in Expression 25 is proportional to the 
density function of a normal random variable Z whose mean is 
u*t and whose standard deviation is a. The last two factors 
are such that their product changes slowly relative to the second 
factor. The consequence of this is that a plot of Expression 25 

against z in units of a for values of t in units of a / u 

has the appearance of the plot of a normal density function for 
values of t greater than four in units of o/u. In particular, 
this indicates that u*t in the exponent of the second factor in 
Expression 25 can be viewed as the radius of an average furthest 
on circle. 

A target's rectangular coordinates and its range and bearing 
from the origin are related by: X = R sin 0 and Y = R cos 0. 

By using these equations, Expression 24 can be transformed to: 

(26) (1/27T) ( r/27TO a ) exp { — [ r a + (ut) a ]/2a a ) Ig (rut/a a ) 

which is the joint density function of the random variables R 

and 0. The marginal density function of R can be obtained 

from Expression 26 by integrating this joint density function 
over the possible values of 0 which in radians is over the 
interval 0 to 2 tt. Consequently, the marginal density for R 
can be obtained from Expression 26 by multiplying it by 2 tt. 
Tabulated values of the cumulative distribution function Fj^(r;t) 
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for the marginal distribution of R for the second case are 
listed in Reference 17. 

Position Distributions Based on the Second Motion Model: In 

the second motion model a target's motion is not restricted to a 
constant course and speed throughout the motion. At the end of a 
time step, a target's course or speed or both can change. In 
general, this leads to a position distribution that can not be 
described analytically. And, in general, a monte carlo 
simulation method is required in order to determine a position 
distribution that is based on such a model. As an example of a 
case in which the distribution can be described analytically, see 
Reference 18. As an example of a case in which it can not, 
suppose a target's initial position is described in terms of a 
number assigned to a subregion in the xy-plane where the number 
assigned represents the probability that the subregion contains 
the target at an initial time. In addition, suppose for each 
subregion a course and speed distribution is determined by 
assigning numbers to course and speed pairs where a number 
represents the probability the target will have the course and 
speed at the initial time given it is in the subregion at that 
time. Next suppose for each course an speed pair there is a time 
distribution that determines the duration of the course and speed 
pair and that the time distribution is determined by a number 
assigned to each discrete time point where the number represents 
the probability that the targets course and speed pair will be 
determined by a new course and speed distribution given it has 



53 



not been determined at an earlier time point. By extending this 
kind of procedure and then implementing it in a monte carlo 
simulation, one can generate complex position distributions that 
describe a target's position at discrete time points in terms of 
probabilities assigned to subregions of a region in which the 
target. 
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XIV. Position Distributions That Change with Search 

Suppose a target's position at some time is described 
by a position distribution. Now suppose information becomes 
available that a search has been conducted for the target and 
that the target has not been detected. Or suppose the 
information is that the target has been detected. In the first 
case, negative information is available that can be used to 
modify the position distribution. In the second case, positive 
information is available that can be used to modify the position 
distribution. In both of the cases, a target's position 
distribution is assumed to be specified in terms of a set of 
discrete probabilities where each probability corresponds to a 
subregion of the region that contains the target and each is the 
probability that the target is in that subregion. 

Position Distributions and Negative Information: For a 

region that contains a target and consists of n subregions, let 
the event = {the target is in the i th subregion}. And let 

the event C = {no detection during a search of the region). 

Then, given an unsuccessful search in the region, the targets 's 
position distribution can be modified as follows: 



(26) 


P(Si|c) = P (C | 


S i )P(S i )/P(C) 




where 


i 1 , 2 , ••• , n 


and P(C) = £ P(C|Sj)P(Sj) 


with the sum 


from 


j = 1 , 2 , ••• , n . 


Note that Equation 26 can 


be obtained by 



using Bayes theorem. To illustrate how Equation 26 might be 
used, suppose that a search in a subregion is considered to be a 
random search and that the sweep width of detection system 
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against the target depends on subregion being searched. For this 
case, let be the area of the i th subregion and let be 

the sweep width in that subregion. Then, given an unsuccessful 
search, P(C|Si) = exp (-Wj_ • li/AjJ where 1^ is the track length 
of the searcher in the i th subregion. Given values for P(S^) , 
$i and 1^ for i = 1,2, ••• ,n a position distribution can be 

determined that has been modified by the negative information. 

Position Distributions and Positive Information: As above, 

let Si the event that a target is in the i th subregion of a 
region that contains the target. However, in the case of 
positive information, instead of the event C occurring, the 
event C = {a contact} occurs. This event is the union of two 
mutually exclusive and exhaustive events: T = {a true contact} 

and F = {a false contact}. Relative to the Venn diagram of 
Figure 1, for a subregion i, Si corresponds to H^, C 
corresponds to Di , T corresponds to (Dj n Hi) and F 
corresponds to (Dj n Hg) . For generality, suppose true contacts 
do not localize a target to a single subregion. Then after a 
search of a region that has resulted in a contact, the target's 
position distribution can be modified as follows: 

(27) P(Si|C) = P (Si | T) P(T|C) + P(Si|F) P(F|C) 

where i = 1,2, ••• ,n. Equation 27 results since T = T n C, 

F = F n C and P(Si|C) = (P[(Si n T) + P(Si n F) ] }/P(C) . The 
probability p = P(t|C) has been called the credibility of the 
contact. In terms of p Equation 27 becomes: 

(28) P(Si|T) = P(Si|T)p + P(Si|F)(l - p) . 
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For the model, P(S-jjT) can be determined by the coverage 
characteristics of the detection system that is used to make the 
contact. In particular, P(SjjT) = P(T| S^) P(S^)/P(T) where 
P(T) = ZP(T| Sj ) P(Sj ) and the sum index j = 1,2, ,n. Then, 

using the correspondence between the events and and the 

events T and ( n ) one can write (Pd) i = p ( T l s i) where 
(Pd)i is an average probability of detection over the i th 
subregion. In one positive information model, P(S-jjF) = P(Sj^), 
the probability the target is in the i th subregion before the 
contact information was available, and p is determined 
subjectively based on factors associated with a search. The 
probability determined by 1 - p has been called the false alarm 
probability. However, it is not the probability of a false alarm 
as it is defined in Section I. Since 1 - p = P(F|c), it 
corresponds to P(H 0 |D 1 ) not to pf = P(DjJh 0 ). 
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XV. Search Models and Search Theory 

Search theory provides a basis for determining optimal 
search plans for a target whose motion and location are 
determined within some bounds. Here, an optimal search plan is 
one for which the probability of finding a target within a given 
length of time is a maximum, the expected time to find a target 
is a minimum given the target is found or a search plan for which 
some other optimal search criterion is satisfied. 

Search theory results are based on models of the search 
process. To the degree that a search model describes a search 
process, an optimal search plan for a target that is based on the 
search model should provide guidance for the development of an 
operationally feasible search plan. However, because of the 
limitations of analytical search models, an optimal search plan 
that is based on an analytical search model may give only initial 
guidance in this regard. The optimal search plans that are 
described below illustrate this. The search plans are based on 
the random search model. Because of this, the requirement on the 
location of search track segments is not realizable and the time 
to resolve false alarms is ignored. 

Optimal search plans based on search models implemented 
through a monte carlo simulation are not considered here. 

However, with sufficient information, such plans have the 
potential of being both implementable and more optimal in a real 
sense than an optimal search plan based on an analytical search 
model . 
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Three Optimal Search Plans: The three optimal search plans 

differ through their definition of optimality. However, each one 
is based on the following search model: A target is fixed at 
some point in a region that consists of n subregions. A search 
in a subregion is a random search in the sense of the definition 
in Section X and a searchers sweep width there is a constant. 

In addition, a search of a subregion will not detect a target 
which is in another subregion. To determine a plan, let 
S^ = {the target is in subregion i) for i = 1,2, ••• ,n and 
let pi = P(Si) be the prior probability that the target is in 
the i*-h subregion. Let W-^ be the sweep width in the i*-* 1 
subregion. Let = Aj/W^ where A^ is the area of the i^h 

subregion and S ^ is the expected track length to find the 
target by a search of the i th subregion given the target is in 
the i^h subregion, a characteristic length. The probability P 
that the target will be detected by a random search is given by: 



where the sum index i = 1,2, ••• ,n and 1 ^ is the track 



The first optimality criterion is: Choose lj_ so that P 

is a maximum subject to the two constraints: 1. 1=21^ and 
2. 1^ > 0 where the index i = 1,2, ••• ,n. Determining this 
choice is a nonlinear optimization problem whose solution is 
given in Reference 19. It is: 



(29) 



P = S [1 - exp (-li/6 i) ] - Pi 



length of the search in the i^h subregion. 



(30) 



li/Si = ln(p i /6 i ) - L ( k) 

lf/Si = 0 



i 1,2, * * * , k 



i = k+ 1 , k+2 , 



n 
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where L(k) = (1/ESj ) j *ln(pj/6j ) ] + 1/26 j and the sum index 
j = 1,2, ,k, where the subregions are relabeled so that the 

following order relation holds: Pi/^i > P2/ 6 2 > ••• > p n /6 n and 
where k is chosen so that for k+1 the solution for l]<+i 
using L(k+1) is either negative or zero. 

The second optimality criterion is: Choose lj^ so that P 

is a maximum subject to the two constraints: 1. c = 2 c-^ and 

2. c^ > 0 where the index i = 1,2, ••• ,n, c^ = k-pl^ is the 
cost of the search in the i^h subregion and k^ is the cost per 
unit track length in that subregion. For this criterion, the 
solution to the corresponding nonlinear optimization problem can 
be obtained from Equation 30 by replacing 6-^ by = k^*&j_ 

and labeling the subregions so that Pi/^i > P2/ e i > *** > Pn/ e i* 
The basis for this can be seen by replacing li/^i by its 
equivalent Cj./ej. in the exponential term in Equation 29. 

The third optimality criterion is: Choose 1^ so that the 

expected utility of finding the target is a maximum subject to 
the two constraints: 1. 1 = Z lj_ and 2. 1^ > 0 where the 
index i = 1,2, ,n. For this criterion, the solution to the 

corresponding nonlinear optimization problem can be obtained from 
Equation 30 by replacing p-^ by q^ where q^ = u^*Pi and u^ 
is the utility of finding the target given it is in the i*-* 1 
subregion. And, in addition, labeling the subregions so that 
qi/6i > q2/£>2 > *** > ^n/^n* The basis for this can be seen by 
multiplying the summation term in Equation 29 by u^ so that 
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Equation 29 gives the expected utility of finding the target 
with the search- 

Equation 30 can be used to determine an order of search for 
the subregions which will effectively minimize the expected track 
length required to detect a target given it is detected. To do 
this, divide the available track length 1 into units small 
enough so that with a single unit only the 1 st subregion would be 
searched. Then allocate one unit to the search of the 1 st 
subregion. If the search is unsuccessful, determine the optimum 
allocation for two units. Then search with a second unit so that 
the first search with the first unit plus the second search with 
the second unit satisfy the optimum allocation for two units. If 
the search is unsuccessful, continue in this fashion until either 
the target is found or all the track length is expended. That 
this allocation order will effectively minimize the expected 
track length required to detect a target given it is detected can 
be argued as follows: Let L be the track length at detection, 

let l u be a unit of track length and let n be the number of 
units. Then the value of the probability P(L < i*l u ) that the 
target will be detected on or before the i th step of the search 
for the given allocation order will be greater than or equal to 
its value for any other allocation order with the same allocation 
step size. Since the value of P(L < 1) will be equal to its 
value for any other allocation order of the optimum allocation 
and since P(L < i*l u |L < 1) = P(L < i*l u )/P(L < 1), the value of 
the distribution function FL(i*l u |L < 1) = P(L < i*l u |L < 1) 
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will also be greater than or equal to its value for any other 
allocation order. This implies that the expected track length 
given detection E(L|L < 1) = £ [1 - FL(i*l u |L < 1) ] where the 
sum index i = 1,2, ••• ,n is effectively a minimum for the 

given allocation order. A search based on the optimum allocation 
given by Equation 30 and the given allocation order is 
equivalent to the following search prcedure: After an allocation 

of track length l u and an unsuccessful search, new values for 
P(S^) are calculated using Equation 26 and then Equation 30 is 
used with these new values to determine the next optimum 
allocation. A discussion of this procedure is given in Reference 
6. And an example of its application is given in Reference 20. 

Equation 30 also defines an optimal search plan for a 
detection system that searches beams and can be described by 
Equation 28 by replacing 1^ by t^ where t^ is the time the 
ith beam is searched and by replacing S ^ by 7^ where t-^ is 
the expected time to detect the target by a search of the i*-* 1 
beam given the target is in the i^-* 1 beam, a characteristic time. 

For a more extensive discussion of search theory and its 
application to military operations research, see Reference 21. 
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